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Abstract: We report on the discovery of high-energy (HE; E > 0. 1 GeV) and very high-energy (VHE; E > 100 GeV) 
7-ray emission from the blazar RBS 0413 with Fermi Large Area Telescope (LAT) and VERITAS. The spectral energy 
distribution (SED), including contemporaneous X-ray and optical observations, is presented. Synchrotron self-Compton, 
external-Compton and lepto-hadronic models are applied to the SED and the results are discussed. 

Keywords: Observations, active galactic nuclei, blazars. 



1 Introduction 

Blazars are active galactic nuclei (AGN) that have their jet 
axes oriented at a small angle with respect to the observer. 
They are broadly classified as either flat-spectrum radio 
quasars (FSRQ) orBL Lacertae objects (BL Lac) according 
to their optical spectra, and are known to emit non-thermal 
radiation characterized by a double-peaked SED. The low- 
energy peak, covering the optical to UV/X-ray bands, is 
usually explained as due to synchrotron emission from rel- 
ativistic electrons in the blazar jets. The origin of the high- 
energy peak, occurring in the X-ray to 7-ray regime, is still 
not completely resolved and could be due to leptonic Q] or 
hadronic |2| processes. RBS 0413 was discovered in the 
X-ray band (IE 0317.0+1834) during the Einstein Medium 
Sensitivity Survey and was optically identified as a BL 
Lac Q. The object was also detected as a radio emitter 
by the Very Large Array of the National Radio Astronomy 
Observatory g). It has a redshift of 0.190 GJ 0. Hav- 
ing a "featureless" optical spectrum [6| and an estimated 
synchrotron peak frequency log(i/ pca k/lHz) = 16.99 [7|, 
RBS 0413 is classified as a high-frequency-peaked BL Lac 
object (HBL; ^ pcak > 10 15 Hz, ®). The MAGIC Collab- 
oration observed the object in 2004 - 2005 and reported a 
7-ray flux upper limit of 4.2 x 10~ 12 ergs cm~ 2 s _1 , at 
200GeV, assuming a T = —3.0 power-law spectrum (9). 
VERITAS observed the source in 2008-2009 season and 
obtained a marginal significance of ~ 3 standard devia- 
tions (cr). In 2009, Fermi-LAT reported HE emission from 
the direction of RBS 0413 flU, triggering new VERITAS 



observations. In October 2009, VERITAS discovered 7-ray 
emission from RBS 0413 ifTTI 



2 Data Analysis and Results 
2.1 VERITAS 

VERITAS is a ground-based VHE 7-ray observatory con- 
sisting of four atmospheric-Cherenkov telescopes, located 
at the Fred Lawrence Whipple Observatory (FLWO) near 
Amado, AZ, USA. It is sensitive to 7 rays with energies 
from 100 GeV up to 30 TeV, with an energy resolution 
of 15-25%. In summer 2009, one of the telescopes was 
relocated and a new mirror-alignment system was intro- 
duced lfl2l . improving the sensitivity of the array IT3l . 
VERITAS can detect a 7-ray source with a flux of 1% of 
the Crab Nebula flux in less than ~ 30 hours, with an- 
gular resolution less than 0.1° and pointing accuracy bet- 
ter than 50". Data analysis consists of a number of steps, 
including calibration, image parametrization, event recon- 
struction, background rejection and signal extraction, as 
described in [14]. 

Twenty-six hours of good-quality VERITAS data on 
RBS 0413 taken between Sep 2008 and Jan 2010 resulted 
in an excess of 108 candidate 7 rays and a 5.50 a detec- 
tion of the blazar. The spectrum above 250 GeV can be 
described with a power law: dN/dE — F (E/E )~ a , 
where F Q = (1.38 ± 0.52) x 10~ 7 TeV" 1 m~ 2 s -1 , 
a = 3.18±0.68 stQt and£o = 300 GeV, with a value of the 
X 2 per degree of freedom (x 2 /dof) of 0.14/2, correspond- 
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Figure 1: RBS 0413 VERITAS (top) and Fermi-LAT (bot- 
tom) light curves. The dashed lines show the average flux. 
The shaded areas in the Fermi-LAT light curve represent 
the time intervals excluded due to solar contamination. 

ing to a fit probability of 93%. No evidence for variability 
is detected (see Figure 1 top panel). 

2.2 Fermi-LAT 

Fermi-LAT is a pair-conversion 7-ray detector sensitive to 
photons in the energy range from below 20 MeV to more 
than 300 GeV [ 15 1. The present analysis includes the data 
taken between 4 Aug 2008 and 4 Jan 2011, which covers 
the entire VERITAS observation interval. "Diffuse class" 
events with energy between 300 MeV and 300 GeV from a 
square region of side 20° centered on RBS 0413 and with 
zenith angles < 100° were selected for this analysis. The 
cut at 300 MeV was used to minimize larger systematic er- 
rors at lower energies. The time intervals when the source 
was close to the Sun (MJD 54954-54974 and 55320-55339) 
were excluded. The data were analyzed with the LAT Sci- 
ence Tools version v9r20p(fl and the post-launch instru- 
ment response functions P6_V1 1_DIFFUSE. The binned 
maximum likelihood tools were used for significance and 
flux calculation lfl6l [TTl . Sources from the 1FGL catalog 
ifTOl located within a square region of side 24° centered 
on RBS 0413 were included in the model of the region. 
The background model includes the standard galactic and 
isotropic diffuse emission component^. 

A point source positionally consistent with RBS 0413 is 
detected with a significance of more than 9 a (test statis- 
tic, TS = 89). The photon energy spectrum is well de- 
scribed by a power law. Replacing the power-law model 
with a log-parabola model does not significantly improve 
the likelihood fit (the likelihood ratio test rejects the power- 
law model with a significance of only 0.35 a). The time- 
averaged integral flux is / (300 MeV < E < 300 GeV) = 



(1.64 ±0.43, 



0.31, 



0.22 s „ fl ) xlO" 



presented in 1101 . In the energy range 100-300 GeV, no 
detection was obtained (TS < 9) and an upper limit at the 
95% confidence level was derived. No evidence for vari- 
ability was found when the data were analyzed using ~ 6- 
month wide time bins (see Figure 2 bottom panel). 

2.3 Swift 

The VERITAS detection triggered a Swift fljD target-of- 
opportunity observation of RBS 0413 on 11 Nov 2009, 
with a total exposure of 2.4 ks. The Swift-XRT photon 
spectrum is described by an absorbed power-law model 
with a fixed galactic hydrogen column density Nr = 
8.91 x 10 20 cm" 2 (19), yielding a x 2 /dof value of 25.9/26. 
Over the energy range 0.3-10 keV, the best-fit photon in- 
dex is T = 2.22 ± 0.07, and the 1 keV normalization is 
(3.31 ± 0.22) xlO keV" 1 m~ 2 s -1 . No flux variability is 
evident over the 2.4 ks exposure. 

UVOT observations were taken in the photometric band 
UVM2 (2246 A) J20). The uvotsource tool was used to 
extract counts, correct for coincidence losses, apply back- 
ground subtraction, and calculate the source flux. The stan- 
dard 5" radius source aperture was used, with a 20" back- 
ground region. The source flux was dereddened using the 
interstellar extinction curve in ETI . The measured flux is 
(2.75 ±0.11) xl0- 12 ergcm" 2 s" 1 . 

2.4 MDM 

The R-band optical data were taken with the 1.3 m 
McGraw-Hill telescope at the MDM observatory on Kitt 
Peak, Arizona, between 10 and 13 Dec 2009. All frames 
were bias-corrected and fiat-fielded using standard routines 
in IRAF [22], and instrumental magnitudes of RBS 0413 
and six comparison stars in the same field of view were ex- 
tracted using DAOPHOT J23) within IRAF. Physical mag- 
nitudes were computed using the physical R-band magni- 
tudes of the six comparison stars from the NOMAD cata- 
log 11241 . assuming that the magnitudes quoted in that cata- 
log are exact. The data show up to ~ 30% daily flux vari- 
ability, with an average magnitude of 17.27. The magni- 
tudes were then corrected for galactic extinction using ex- 
tinction coefficients calculated following 11251 . taken from 
NEE0 and then converted into v¥ v fluxes. 

3 SED Modeling 

The non-thermal continuum of RBS 0413 exhibits a 
double-peaked shape, typical for blazars. In this study, 
three different time-independent models were used to de- 
scribe the observed SED, using contemporaneous X-ray, 
UV and optical (R-band) data complementing the Fermi- 
LAT and VERITAS observations (see Figure 2). All model 



The spectral index is 1.57 ± 0.12 s t at +0.11 
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The spectral points were calculated using the procedure 



1. http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview.html 

2. http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 

3. http://ned.ipac.caltech.edu/ 
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Figure 2: RBS 0413 spectral energy distribution. 



flux calculations were corrected for extragalactic back- 
ground light (EBL) absorption using the model of [ 26] . The 
first modeling, carried out assuming a pure SSC process, 
requires a magnetic field energy density which is only 6% 
of that required for equipartition with the relativistic elec- 
tron distribution. In addition, the optical (R-band) spec- 
trum is not reproduced, the model Fermi spectrum is too 
hard and the model VERITAS spectrum is too soft, albeit 
within the errors in both cases. Next, a combined SSC 
+ EC model was adopted. The addition of an EC com- 
ponent gives a better description of the optical and GeV 
parts of the SED with parameters very close to equiparti- 
tion between the magnetic field and relativistic electrons, 
but slightly underproduces the TeV data. Finally, a com- 
bined lepto-hadronic jet model as described in [27] was 
used. This model gives the best description for the data, 
with the system close to equipartition. On the other hand, 
it is the least constraining model given the large number of 
free parameters. 

4 Summary 

VERITAS and Fermi-LKT detected VHE and HE 7-ray 
emission from the HBL RBS 0413. No evidence for vari- 
ability was seen in either data set. The contemporaneous 
SED of RBS 0413 covering the energy range from the op- 
tical band to TeV 7 rays was studied and three different 
emission models were applied to the data. A pure SSC 
model fails to reproduce the R-band data and requires a 



magnetic field far from equipartition. The addition of an 
EC component yields an improved result but slightly un- 
derproduces the TeV data. The lepto-hadronic model gives 
the best description of the data, with the system close to 
equipartition, but it is the least constraining model since it 
has the largest number of free parameters. 
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